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ABSTRACT 

Water, methane and carbon-monoxide are expected to be among the most abundant molecules besides 
molecular hydrogen in the hot atmosphere of close-in extrasolar giant planets. Atmospheric models for 
these planets predict that the strongest spectrophotometric features of those molecules are located at 
wavelengths ranging from 1 to 10 /im making this region of particular interest. Consequently, transit 
observations in the mid-IR allow the atmospheric content of transiting planets to be determined. We 
present new primary transit observations of the hot-jupiter HD 189733b, obtained simultaneously at 
4.5 and 8 /xm with the Infrared Array Camera (IRAC) onboard the Spitzer Space Telescope. Together 
with a new refined analysis of previous observations at 3.6 and 5.8 /xm using the same instrument, 
we are able to derive the system parameters, including planet-to-star radius ratio, impact parameter, 
scale of the system, and central time of the transit from fits of the transit light curves at these four 
wavelengths. We measure the four planet-to-star radius ratios, to be {Rp/Ri,)^ e/^m = 0.1545 ±0.0003, 
(i?p/i?,)4.5,.m = 0.1557±0.0003, (i?p/i?*)5.8Mm = 0.1547±0.0005, and (i?p/i?08,.m = 0.1544±0.0004. 
The high accuracy of the planet radii measurement allows the search for atmospheric molecular ab- 
sorbers. Contrary to a previous analysis of the same dataset, our study is robust against systematics 
and reveals that water vapor absorption at 5.8 fim is not detected in this photometric dataset. Fur- 
thermore, in the band centered around 4.5 /im we find a hint of excess absorption with an apparent 
planetary radius ARp/R^, =0.00128±0. 00056 larger {2.3a) than the one measured simultaneously at 
8 /im. This value is 4tT above what would be expected for an atmosphere where water vapor is the 
only absorbing species in the near infrared. This shows that an additional species absorbing around 
4.5 /im could be present in the atmosphere. Carbon monoxide (CO) being a strong absorber at this 
wavelength is a possible candidate and this may suggest a large CO/H2O ratio between 5 and 60. 
Subject headings: eclipses planetary systems stars: individual: (HD189733b) techniques: photometric 



1. INTRODUCTION 

The theory of transmission spectroscopy was developed 
in detail by the pioneering works of Seager & Sasselov 
(2000), Brown (2001), Burrows et al. (2001), and Hub- 
bard et al. (2001). This method is now widely used to 
probe the atmospheric structure and composition of the 
transiting planets. Applied to the first transiting extra- 
solar planet known, HD 209458b, this method allowed 
for the first detection of an extrasolar planetary atmo- 
sphere (Charbonneau et al. 2002), and of the discoveries 
of an escaping exosphere (Vidal-Madjar et al. 2003, 2004, 
2008; Ehrenreich et al. 2008), and a hot hydrogen layer 
in the middle atmosphere (Ballester et al. 2007). 

HD 189733b is a transiting hot-jupiter which orbits a 
small and bright main sequence K star, and shows a tran- 
sit occultation depth of ~2.5% (Bouchy et al. 2005). The 
planet has a mass of Mp = 1.13 Jupiter mass (Mj) and a 
radius of Rp = 1.16 Jupiter radius in the visible (Bakos 
et al. 2006b; Winn et al. 2007). The short period of 
the planet (~2. 21858 days) is known with high accu- 
racy (Hebrard & Lecavelier des Etangs 2006; Knutson 
et al. 2009). Sodium in the atmosphere of the planet has 
been detected through ground-based observations (Red- 
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field et al. 2008). Using the Advanced Camera for Sur- 
vey aboard the Hubble Space Telescope (HST), Pont et 
al. (2008) detected atmospheric haze interpreted as Mie 
scattering by small particles (Lecavelier des Etangs et al. 
2008a). 

The red part of the visible to the Infrared (IR) is a 
spectral region of particular interest to detect molecular 
species using transmission spectroscopy (e.g., Tinetti et 
al. 2007). The rotational- vibrational transition bands 
of water (H2O), carbon monoxide (CO), and methane 
(CH4) are anticipated to be the primary sources of non- 
continuum opacity in hot- Jovian planets. Their rela- 
tive abundances are determined by the C/0 ratio (Liang 
et al. 2003, 2004; Kuchner & Seager 2006). Carbon 
monoxide is expected to be the dominant carbon-bearing 
molecule at high temperatures and CH4 the dominant 
one below ~ 1000 K (Kuchner & Seager 2006; Cooper 
& Showman 2006; Sharp & Burrows 2007). In the case 
of HD 209458b's atmosphere. Brown et al. (2002) and 
Deming et al. (2005) reported exploratory attempts to 
detect CO during transit from ground based observa- 
tions. In the same atmosphere, the possible presence of 
H2O (Barman et al. 2005), TiO and VO (Desert et al. 
2008; Sing et al. 2008a,b), and the detection of the main 
constituent, H2, through Rayleigh scattering (Lecavelier 
des Etangs 2008b), have been reported using observa- 
tions in the visible from the Space Telescope Imaging 
Spectrograph aboard HST. 

Richardson et al. (2006) obtained the first infrared (IR) 
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transit measurement for this planet using the Spitzer 
Space Telescope (Werner et al. 2004) at 24 ^m, and 
derived the planet radius with high accuracy. 

HD 189733b has been observed at 8 fj,m (Channel 4) 
with the Infrared Array Camera (IRAC; Fazio et al. 
2004) aboard Spitzer. The planet-to-star radius ratio 
at 8 ^m was found to be (i?p/i?*)8/xm = 0.1545 ± 0.0002 
(Knutson et al. 2007a). Agol et al. (2008) also observed 
this planet at 8 /im during 7 primary transits and used 
the high photometrical precision at this wavelength for 
the searching companions using Transit Timing Varia- 
tions (TTV) method. 

More recently, secondary eclipse measurements re- 
vealed the presence of water, carbon monoxide and car- 
bon dioxide in the dayside emission spectrum of the 
planet, between 1.5 and 2.5 /im (Swain et al. 2009), and 
water absorption with the possible presence of carbon 
monoxide as a contributor near the Spitzer /IRAC 4.5 /im 
(Grillmair et al. 2008). 

Results at 3.6 /im and 5.8 /im (Channel 1 and 3) 
from primary transit observations of HD 189733b using 
Spitzer /IRAC have already been published but led to two 
different conclusions (Ehrenreich et al. 2007; Beaulieu et 
al. 2008). From the analysis of Beaulieu et al. (2008); 
Tinetti et al. (2007) concluded that water vapour is de- 
tected in the atmosphere of the planet, whereas Ehren- 
reich et al. (2007) claim that uncertainties on the mea- 
surements are too large to draw any firm conclusions. In- 
dependently, a detection of water and methane has also 
been obtained using observations between 1.5 and 2.5 /im 
from the Near Infrared Camera and Multi-Object Spec- 
trometer aboard HST (Swain et al. 2008). 

Here we describe Spitzer observations collected during 
two primary transits of HD 189733b, with the intent to 
measure its radius at the four IR bandpasses of the IRAC 
instrument. We analyze, in a consistent way, the obser- 
vations from four bands centered at 3.6, 4.5, 5.8 and 8 
/im. We probe the planet atmospheric composition by 
comparing results from these four photometric bands. 
We first describe the observations and data reduction in 
Sect. [2] In Sect.[3l we describe the model and techniques 
used to estimate the physical and orbital parameters of 
HD 189733b. We discuss our results in Sect.[4]in the hght 
of previous analysis and theoretical predictions. We fi- 
nally summarize our work in Sect.O 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Observations 

We obtained Spitzer Guest Observer's time in Cycle 
3 and 4 (PI: A. Vidal-Madjar; program IDs 30590 & 
40732), in the 4 Spitzer /IRAC bandpasses. Our primary 
scientific objectives were to detect the main gaseous con- 
stituents (H2O and CO) of the atmosphere of the hot- 
Jupiter HD 189733b. First observations of the system 
were performed on 2006 October 31 simultaneously at 
3.6 and 5.8 /im (channels 1 and 3) (Ehrenreich et al. 
2007). The second part of the program was completed 
on November 23, 2007 using the 4.5 and 8 /im channels 
(channels 2 and 4) following the same methods. The sys- 
tem was observed using IRAC's stellar mode during 4.5 
hours for each visit, upon which 1.8h was spent in plan- 
etary transit. The observations were split in consecutive 
subexposures at a cadence of 0.4 s for channels 1 and 2 
and 2.0 s for channels 3 and 4. We obtained a total of 



1936 frames for channels 1 and 3, and 1920 frames for 
channels 2 and 4. 

While planning the observation in the Astronomical 
Observing Request (AOR) format, we carefully selected 
a pixel area avoiding dead pixels. We also purposely did 
not dither the pointing in order to keep the source on a 
given pixel of the detector. Doing so minimizes errors 
from imperfect flat-field corrections and thus increase 
the photometric accuracy. This common observational 
strategy matches that of recent Spitzer observations of 
HD 189733 and GJ 436 (Knutson et al. 2007; Deming et 
al. 2007, Gillon et al. 2007; Agol et al. 2008) 

2.2. Data reduction 

We used the Spitzer /IRAC Basic Calibrated Data 
(BCD) frames in the following analysis. These frames 
are produced by the standard IRAC calibration pipeline 
and include corrections for dark current, flat-fielding, de- 
tector nonlinearity, and conversion to flux units. We first 
find the center of the PSF of the star to a precision of 
0.01 pixel using the DAOPHOT-type Photometry Proce- 
dures, CNTRD, from the IDL Astronomy Library which 
computes the centroid of a star using a derivative search. 
We find that the position of the center of the star varied 
by only less than 10% of a pixel during the whole obser- 
vation. We also tested the position of the centroid of the 
star using GCNTRD, which computes the stellar centroid 
by Gaussian fitting. As a final test, we used a weighted- 
position sum of the flux procedure in a 5 x 5 pixel box 
centered on the approximate position of the star. These 
two last methods give very similar results, which are 
in agreement within 0.02% with the positions extracted 
using CNTRD. We used the APER routine, which we cus- 
tomized to perform a weighted aperture photometry by 
summing the weighted background-subtracted flux, on 
each pixel, within an aperture of a given radius (Horne 
1986; Naylor 1998). We used a radius of 5 pixels to 
minimize the contribution of HD 189733B (Bakos et al. 
2006a), the closest star in the field of view. We checked 
that our results remain the same when using a radius 
varying from 4 to 6 pixels. The background level for each 
image was determined with APER by the median value of 
the pixels inside an annulus, centered on the star, with 
an inner and outer radii of 16 and 18 pixels, respectively. 
This constant level is subtracted from each pixel of a 
subexposure to create the background-subtracted image. 
The background errors are dominated by statistical fluc- 
tuations. The PSF, used for weighting, is estimated in 
each channel as the median of the background-subtracted 
fluxes. The estimated error on the weighted integrated 
flux is calculated as the square-root of the photon-noise 
quadratic sum (Horne 1986; Naylor 1998). The four raw 
weighted light curves obtained with this method are plot- 
ted in the top panel of Fig. [TJ After producing a time 
series for each channel, we iteratively select and trim out- 
liers greater than 3(t by comparing the measurements to 
a transit light curve model. Doing so, we remove any 
remaining points affected by transient hot pixels. We 
discarded 24, 26, 46 and 19 exposures in channel 1 to 4 
respectively, which represent 1 to 2% of the total number 
of data points. 

^ |http : // idlastro . gsf c ■ nasa ■ gov/homepage ■ html | 
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3. FITTING THE SPITZER TRANSIT LIGHT CURVE 

3.1. Out- of- Eclipse Baseline 

Because of instrumental effects, the measured stellar 
flux out-of-transit is not constant, but is seen to vary in 
time. To correct for the instrumental effects, we define 
a baseline function of time for each Channel. The base- 
line function is time- and wavelength-dependent and is 
used to normalize a given time series. We find that a 
linear function of time represents well the data in chan- 
nels 1 and 2. On the other hand, in Si:As based de- 
tectors (channels 3 and 4), the effective gain, and thus 
the measured flux in individual pixels, drift non linearly 
over time. This effect, called detector ramp, is well doc- 
umented (Deming et al. 2005, 2006; Harrington et al. 
2007; Knutson et al. 2007). To correct for the ramp and 
other baseline effects, we adopted a non-linear function of 

time FbaseUne = ^0+^1 xt-^^a log(t-io)+^3 log^(i-^o), 
where FbaseUne IS the flux of the central star HD189733 
without planetary transit and to was fixed to a time a 
few minutes before the first observations. As a check, we 
also tested linear, polynomial and exponential baseline 
functions but find that the logarithmic baseline provides 
the best results for channels 3 and 4 (see also Sect. |4]). 
We find that baselines other than logarithmic introduce 
significant systematic errors when determining planetary 
parameters as seen in Fig. [2] We fit the detector correc- 
tion coefficients simultaneously with the transit-related 
parameters, allowing us to take into account how changes 
in the correction coefficients may impact the transit pa- 
rameters. 

3.2. Pixel-phase effect correction 

Telescope jitter and intra-pixel sensitivity variations 
are also responsible for fluctuations seen in the raw light 
curves (Fig.[l]upper panel), most severe in channel 1. We 
find that the pixel-phase effect is negligible in our dataset 
for the three channels 2, 3 and 4. A description of this ef- 
fect is given in the Spitzer /lYiAC data handbook (Reach 
et al. 2006, p. 50; Charbonneau et al. 2005). The method 
reported by Morales-Calderon et al. (2006) and also ap- 
plied by the two previous analyses of the present dataset 
for channel 1 and 3 (Ehrenreich et al. 2007; Beaulieu et 
al. 2008) consists in calculating a pixel-phase dependent 
correction by fitting the light curve variations induced by 
this effect using a function of one single parameter for the 
X and Y target positions. We find that this single order 
function, described with only one parameter, poorly cor- 
rects the light curves for the pixel-phase effect, and thus, 
systematics remains in the residuals (Sect. |4] and Fig. [3]). 
We conclude that this function (Reach et al. 2006) is 
not appropriate to high precision photometry. To bet- 
ter decorrelate our signal from pixel-phase variations, 
we adopted a different correction, based on a quadratic 
function as described by Knutson et al. (2008) where we 
added a cross term (X5), Fcorr = F{\ -\- Ki{x — xq) -\- 
K2{x-xof+K3{y-yQ)+Ki{y-yof+K5{x-Xf)){y~yo)), 
where xq and yo are the integer pixel numbers contain- 
ing the source centroid and F is the measured flux of 
the central star HD 189733. The Ki coefficients are the 
decorrelation factors that have to be derived from flts to 
the the transit lightcurve (TLC). A quadratic function 
of X and Y , with flve parameters, improves the fit com- 
pared to a linear function, from a of 2077 to 1947 for 



n = 1912 degrees of freedom. Introduction of the cross 
term [K^) also slightly improved the fit. We find that 
adding higher-order terms to this equation does not im- 
prove the fit for channel 1. As an additional test, we 
applied this decorrelation procedure to one of the bright 
stars in the field of view, resulting in a time series that 
showed no significant deviations from a constant bright- 
ness. 

The intra-pixel sensitivity is also expected to contam- 
inate the photometry in channel 2. However, we did not 
notice such a strong effect as seen in Channel 1. We con- 
cluded that the central star was localized on a part of the 
array which has a flat photometric response (pixel refer- 
ence: 147.20,198.25). This part of the detector might be 
of particular interest for the warm Spitzer mission. 

3.3. Transit light curve model with limb darkening (LD) 

corrections 

We parameterized the transit light curve with 4 vari- 
ables: the planet-star radius ratio Rp/Ri,, the orbital 
distance to stellar radius ratio a/ Ri,, the impact param- 
eter b, and the time of mid-transit Tc- We used the IDL 
transit routine OCCULTNL developed by Mandel & Agol 
(2002) for the transit light curve model. 

For each channel, we calculated a theoretical limb- 
darkening model (Kurucz 1979) with TeB = 5000K, 
log (7 = 4.5, [Fe/H] solar, and fit this model^ to derive the 
4 non-linear limb-darkening coefficient defined by Claret 
(2000) and presented in Table[TJ We found that account- 
ing for the effects of limb-darkening decreased the re- 
sulting best-fit transit depth by 0.5cr at 8.0 /xm to la at 
3.6 /im. We tested the robustness of our result for several 
limb-darkening corrections. We derived planet-to-star 
radius ratios using linear and quadratic limb darkening 
coefficients. In these cases the resulting best-fit transit 
depth decreased by OAa for channel 1 and 2 and by 0.8a 
for channel 3 and 4 compared to a 4th degree polynomial 
limb-darkening corrections. Finally, we modified the val- 
ues of the non-linear limb-darkening coefficients (Claret 
2000) presented in Table [T] by 10% and found that the 
planet-to-star radius ratios changed by less than O.Str. 

3.4. The fitting procedure 

We performed a least-squares fit to our unbinned data 
over the whole parameter space {Rp/R^,, a/Ri,, b, Tc, 
Ai, Ki). In order to find the best-fit observables to 
the data, we used the MPFIT package^, which performs 
a Levenberg-Marquardt least-squares fit. We combined 
the baseline function and the pixel-phase correction de- 
scribed above with the transit light curve function which 
takes limb darkening correction into account. We ap- 
plied this method on the four channels independently to 
calculate the four observables. 

3.5. Mean values and errors determination 

We used a bootstrap method to determine the mean 
value, the statistical and systematical errors for the mea- 
sured parameters. The possible presence of correlated 
noise in the light curve has to be considered (Pont et al. 
2006). 

^ See|http : //kurucz .harvard. e du/gr Ids / | 

® |http : //cow .physics ■ wise . edu /^crai gm/ldl/idl ■ html | 
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The bootstrap technique we used take into account 
both the red and white-noise. To estimate the system- 
atical errors due to intrapixel sensitivity and basehne 
corrections, we randomly padded the beginning of the 
transit light curve from different phase before -0.03. Few 
thousand transit light curves with different baseline du- 
ration are produced that way. Additionally, we randomly 
varied each photometric measurement, within their esti- 
mated error bars, following a normal distribution in or- 
der to derive the statistical errors on the derived param- 
eters. Totally, 4000 synthetic transit light curves were 
produced. We then fitted (Sect. 13. 4p these transit light 
curves to derive a new set of parameters and to extract 
their means and their corresponding la statistical and 
systematical errors. 

As an additional test for the errors, we measured the 
errors using the Prayer Bead method (Moutou et al. 
2004, Gillon et al. 2007). In this case, the residuals 
of the initial fit are shifted systematically and sequen- 
tially by one exposure, and then added to the eclipse 
model before fitting. The purpose of this procedure is to 
take into account the actual covariant noise level of the 
light curve. Using this method, we found negligible red- 
noise after corrections and obtained errorbars equivalent 
to the systematical and statistical errors derived with the 
bootstrap method presented above. 

We examined the residuals from the best fit of each 
synthetic transit light curve and for each channel inde- 
pendently. We measured a RMS residual between 2.5 
and 3.3 x 10"'^ on normalized flux for all exposures, de- 
pending on the channel (See case of Channel 4 bottom 
panel Fig. The RMS residuals are 20% larger than 
the expected photon-noise, and stays constant over the 
duration of the observation. We find that the scatter of 
residuals in binned exposures decreases with bin size as 
jV~i/2 for bins of up to 150 points for all the channels, as 
expected for photon- noise (See channel 1, bottom curve 
in Fig. [3]). 

4. RESULTS AND DISCUSSION 

The results of the independent fits are given in TableO 
The most interesting parameters are the 4 planet-to-star 
radius ratios which are discussed in this section. 

4.1. Comparisons with previous results 

The main difference between our analysis of channels 1 
and 3 and the previous studies (Ehrenreich et al. 2007, 
Beaulieu et al. 2008) resides in the pixel-phase decorre- 
lation and baseline functions adopted. 

In channel 1, the flux measured is strongly correlated 
with the position of the star on the detector array. We 
tested the influence of the pixel-phase decorrelation func- 
tion on channel 1 by comparing a correction function 
with one parameter, as used in the previous studies 
(Ehrenreich et al. 2007; Beaulieu et al. 2008), to a 
2nd degree function with 5 parameters as described in 
Sect. 13.21 For each corrected photometric time-series, 
the red noise was estimated as described in Gillon et al. 
(2006), by comparing the rms of the unbinned and binned 
residuals. In the case of channel 1, the fit is largely im- 
proved by fitting 5 parameters compared to only one pa- 
rameter for the X and Y star position and only negligi- 
ble systematics effects remain in the residuals (Fig. [3]). 
Contrary to the two previous analysis, removing photo- 



metric points which have an extreme pixel-phase value 
does not change our result, since they are well corrected 
by the present decorrelation methods. We thus kept all 
the photometric measurements. 

The baselines are known to be inherently linear for 
channel 1 and 2 and logarithmic for channel 3 and 4 
(Knutson et al. 2008). However, we tested three different 
baselines functions: polynomial of one, two and three de- 
grees, an exponential with a polynomial, and logarithmic 
with a linear function as describe in Sect. 13. ll We tested 
the robustness of each of these different baselines using 
the same test as described in Ehrenreich et al. (2007), 
i.e. by dropping the first exposures from the beginning 
of the observations (exposure with a phase smaller than 
-0.03). We found that for channels 3 and 4, the best 
fits are obtained when using a logarithmic baseline (see 
lower panel in Fig. In the case of a linear, polyno- 
mial or exponential baseline, the fitted parameters show 
large variations with the number of data point removed 
(Fig. [21 upper panel). The radius dramatically changes 
according to the number of points removed when using a 
linear or a third degree polynomial baseline functions in- 
dicating that systematics errors remain in the corrected 
data. In the case of a logarithmic baseline, the radius 
extracted does not depend on the number of removed 
points. The logarithmic function is the only one which 
allows the observable parameters to oscillate around the 
same value independently of removed exposures. Only 
small systematics still remain and they are included in 
the final error bar. Consequently, our estimates of the 
error bars on the measured parameters in channel 3 and 
4 are conservative; they are larger by about 30% than 
the error bars due to the photon noise. 

The present study shows that Ehrenreich et al. (2007) 
and Beaulieu et al. (2008) did not apply a sufficient pixel- 
phase decorrelation for channel 1. Consequently, sys- 
tematics remained in these two previous studies (Fig. [3]). 
Furthermore, none of these studies did use a logarith- 
mic baseline for channel 3. This explains why Ehrenre- 
ich et al. (2007) obtained larger systematics error bars 
in this channel 3 when exploring this issue. Our analy- 
sis also shows that by using linear baselines, the result 
of Beauheu et al. (2008) at 5.8 /zm is affected by large 
systematics which are not included in their error bars es- 
timations. Having understood the reasons of the discrep- 
ancies, the results in the present paper can be considered 
with most confidence. 

The planet-to-star radius ratio at 8 /im that we derive 
(Fig.|4]and Table[2]) is in agreement with a previous mea- 
surement obtained with the same instrument and same 
channel, but in subarray mode (Knutson et al. 2007a). 
The error bars we obtain are slightly larger, since they 
also includes systematic errors from the baseline (Fig. [2]). 

4.2. Discussion 

We derived the fitted parameters and their error bars 
in a consistent way. From our fits, we evaluated at each 
wavelength the radius ratio i?p/i?*, the impact parame- 
ter b, the system scale a/i?^ and the central time of the 
transit (Table [2]). 

The measured central time of the transit can 
be compared to the expected transit time from 
known ephemeris. Several observations of transits of 
HD 189733b have been published recently with improved 
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ephemeris (Bakos et al. 2006b; Winn et al. 2007a; Knut- 
son et al. 2007; Pont et al. 2007; Knutson et al. 2009). 
We find that the central times that we obtained are in 
agreement within the error bars for each set of simulta- 
neous observation, between channels 1 and 3 for the first 
epoch and channels 2 and 4 for the second one. 

Although data from different channels have been fitted 
independently, our measurements of the impact param- 
eter and of the system scale both show the same be- 
havior as a function of the channel: they are consistent 
between channels observed simultaneously, but disagree 
between channels observed at the two different epochs. 
This strongly suggests that unrecognized systematics re- 
mains between the observations obtained at two epochs; 
these systematics could be of instrumental or more likely 
of astrophysical origin (e.g, star spots). We therefore 
consider that comparison of radius ratios should pref- 
erentially be made between measurements taken at the 
same epochs. 

The impact parameter and the system scale values 
can also be compared with the most accurate measure- 
ments (Pont et al. 2007; b = 0.671 ± 0.008 and a/R^ = 
8.92 ±0.009), where the degeneracy between primary ra- 
dius and orbital inclination could be lifted. Our results 
from our primary observations at 3.6 /im and 5.8 /im are 
in agreement with previous analysis of the same data set 
(Ehrenreich et al. 2007; Beaulieu et al. 2008), but are in- 
consistent at 3tT level with the Pont et al. (2007) values. 
This suggests that the measurements in these two chan- 
nels may be affected by, for instance, stellar spots. The 
measurements of the impact parameter and the system 
scale in channels 2 and 4 data at 4.5 /im and 8 fim are in 
agreement with the measurements by Pont et al. (2007). 
This reinforces the confidence in our determination of the 
physical parameters for channels 2 and 4. 

The planet-to-star radius ratio measured at 3.6 fxrn 
is larger than the value extrapolated from simultane- 
ous measurements at 5.8 fim assuming that only water 
molecules contribute to the absorption in this channel. 
Therefore, other species must be present in the atmo- 
sphere of the planet and absorb at this wavelength. The 
possible presence of methane (CH4) could also contribute 
to the opacity at 3.6 ^m. Interestingly, the planetary ra- 
dius measured at this wavelength is in agreement with 
the extrapolation of the measurements obtained with 
HST/ACS in the visible (Pont et al. 2007, 2008) and 
assuming Rayleigh scattering as needed to interpret the 
large variations of the planet radius between 0.55 and 
1.05 /xm (Lecavelier des Etangs et al. 2008b). 

Comparing the planet-to-star radius ratio obtained si- 
multaneously in channels 2 and 4, we find a difference 
(Ai?p/i?^)4.5-8 /.m = 0.00128 ± 0.00056. This corre- 
sponds to a 2.3cr signature of an excess absorption at 
4.5 fim relative to 8 fiui. The radius at 4.5 /im is also 
4(7 above the expected value if only water molecules 
were contributing to the absorption in this channel (see 
Sect. 14. 3p . At this wavelength, another absorbent is 
needed. The new observations at 4.5 /xm are not af- 
fected by the systematics one faces in at 3.6 /xm, mainly 
the pixel phase effect, since the central star was localized 
on a part of the array which has a flat photometric re- 
sponse. Hence data reduction and analysis of data taken 
in channel 2 at 4.5 fim are straight forward compared 
to data in the 3 other channels. As a consequence, the 



parameter extraction in this channel is the most robust 
of the 4 channels; conclusions drawn from this channel 
can be considered with more confidence. 

4.3. Possible CO signature and large C/0 ratio scenario 

From Fig.lHand Table[21 a possible excess of absorption 
could be present in channel 2 around 4.5 fim by compar- 
ison with channel 4. At this wavelength, all main atmo- 
spheric constituents (H2, H2O, CH4) have no strong spec- 
tral features, except CO (Sharp & Burrows 2007; Fig 3). 
We therefore conclude that CO could be a candidate for 
identification of this possible absorption signature. 

Following Eq. 2 of Lecavelier des Etangs et al. (2008a), 
if ARp is the variation of the apparent planetary radius 
between channels 2 and 4, we have 

ARp=H\n^, 

where H is the atmospheric scale height, ^ and a are 
the abundance and cross-section of the main absorbent 
in channels 2 and 4, respectively. For a tempera- 
ture T, H is given hy H — kT/fig, where /x is the 
mean mass of atmospheric particles taken to be 2.3 
times the mass of the proton, and g the gravity. For 
HD 189733b, we have i?/i?*=0.00045 at 1500 K. With 
Ai?p/i?*=0.00128±0.00056, we find 



logio^ = 1.245 ± 0.540 ' 



T 



^4(74 



V1500K 



The H2O absorption cross-section is smaller in chan- 
nel 2, o-H2o(4.5 /xm) = 2 x 10"^^ cm^/molecule, than 
in channel 4, (7h2o(8 /xm) — 2 x 10^^° cm^ /molecule 
(Sharp & Burrows 2007). Therefore, assuming that H2O 
is the main absorber at 8 /xm, we find that the planet-to- 
star radius ratio at 4.5 /xm should be (i?p/i?*)4.5 — 
0.1534±0.0005. This value is 4cr below our effective mea- 
surement. Therefore, species other than H2O should be 
present. Because CO has a strong absorption signature 
at this wavelength, it is a possible candidate for explain- 
ing the larger planetary radius measured at 4.5 /xm. 

CO has a similar cross section at 4.5 /xm as wa- 
ter at 8 /xm : crH2o(8/xm) = crco(4.5/im) = 2 x 
10^^° cm^/molecule, therefore, the two measurements in 
channels 2 and 4 lead to the conclusion that 



logio^co/CH20 = 1.245 ±0.540 



/ T 



V1500K 



This corresponds to CO abundance of 5 to 60 times 
the H2O abundance. Note that the lower limit of the 
CO/H2O ratio 6) is about 3.5(J higher than the high- 
est ratio at which CO absorption at 4.5 /im are overcome 
by H2O absorption (CO/H2O = 0.1). In other words, at 
CO/H2O ratio above = 0.1, CO is the main absorbent 
and can be detected at 4.5 /xm; our measurements show 
an absorption 3.5(7 above this level. 

AC/0 ratio above the solar value could enhance, at 
high temperature, the CO abundance with respect to 
the abundance of II2O in the atmosphere of HD189733b 
(Tinetti et al. 2007). Indeed, carbon-rich planetary 
environments do exist: a C/0 ratio much larger than 
solar has been measured with FUSE and HST in 
the planetary system of f3 Pictoris (Jolly et al. 1998; 
Lecavelier des Etangs et al. 2001; Roberge et al. 
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2006). Planets evolving in such environments could 
have a high C/0 ratio in their atmospheres. The large 
CO/H2O ratio would appear at odd with the detection 
of CH4 claimed by Swain et al. (2008), though the 
presence of CO is more consistent with the temperature 
above 1000 K as estimated from primary and secondary 
transit observations. Furthermore, simulations show 
that vigorous dynamics caused by uneven heating of 
tidally locked planets can homogenize the CO and CH4 
concentrations (Cooper & Showman 2006). However, 
recent 3D models (Showman et al. 2008) shows that the 
planet could experience dramatic day/night temperature 
variations along the planetary limb. As a consequence, 
the detection of both CO and CH4 molecules by primary 
transit observations could be compatible if they are 
located in different limb regions. 

5. SUMMARY 

We have performed and analyzed new Spitzer /IRAC 
primary eclipse observations of HD 189733b at 4.5 fim 
and 8 /xm and reanalyzed observations at 3.6 //m and 
5.8 ^m in a consistent way. Our analysis is more robust 
against systematics than the previous ones (Ehrcnrcich 
et al. 2007, Beauheu et al. 2008), especially regarding 
the pixel-phase effect and baseline correction. Systematic 
error corrections can now be applied in all channels. 

As a consequence, contrary to what was claimed in a 
previous study, we do not find the excess absorption at 
5.8 ^m compared to 3.6 fim that was initially found by 
Beaulieu et al. (2008) and interpreted by water absorp- 
tion by Tinetti et al. (2007). Therefore, other species 
such as methane must be present in the atmosphere of 
the planet and absorb at this wavelength. Noteworthy, 
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the measured radius at 3.6 /im is compatible with the ra- 
dius extrapolated assuming Rayleigh scattering absorp- 
tion by small particles. 

Interestingly in channel 2, the star is centered in a 
region of the array which has a flat photometric response. 
Therefore, the resulting measurements in this channel are 
robust. 

Comparing the planet-to-star radius ratios obtained si- 
multaneously at 4.5 and 8 /zm, we noticed a 2.3a excess 
absorption at 4.5 /im. Furthermore, the radius ratio at 
4.5 /im stands 4ct above the expected value if only wa- 
ter molecules were contributing to the absorption. If this 
signature is confirmed, it could be interpreted by absorp- 
tion by CO molecules. In that case, we estimate a ratio 
of CO/II2O abundances of 5 to 60, possibly indicating a 
high C/0 ratio. 

The Near Infrared Spectrograph (NIRSpec) aboard the 
future James Webb Space Telescope will enable medium- 
resolution spectra over a wavelength range from 1 to 
5 /im and thus allow for clear identification simultane- 
ously and unambiguously of H2O, CII4 CO and CO2. 



We thank the referee who greatly contributed to im- 
prove the paper, as well as H. Knutson and Ph. Nutz- 
man for helpful discussion on data reduction and anal- 
ysis. D.K.S. is supported by CNES. This work is based 
on observations made with the Spitzer Space Telescope, 
which is operated by the Jet Propulsion Laboratory, Cal- 
ifornia Institute of Technology under a contract with 
NASA. D.E. acknowledges financial support from the 
French Agence Nationale pour la Recherche through 
ANR project NT05-4.44463. 



Desert, J.-M., Vidal-Madjar, A., Lecavelier Des Etangs, A., Sing, 

D., Ehrenreich, D., et al. 2008, A&A, 492, 585 
Ehrenreich, D., Lecavelier des Etangs, A., Hebrard, G., et al. 2008, 

A&A, 483, 933 

Ehrenreich, D., Hebrard, G., Lecavelier des Etangs, A., et al. 2007, 
ApJ, 668, L179 

Ehrenreich, D., Tinetti, G., Lecavelier des Etangs, A., Vidal- 
Madjar, A., & Selsis, F. 2006, A&A, 448, 379 
Hc5brard, G., & Lecaveher des Etangs, A. 2006, A&A, 445, 341 
Fazio, G. G., et al. 2004, ApJS, 154, 10 
Fortney, J. J., & Marley, M. S. 2007, ApJ, 666, L45 
Gillon, M., Pont, F., Demory, B.-O., et al. 2007, A&A, 472, L13 
Gillon, M., Demory, B.-O., Barman, T., et al. 2007, A&A, 471, L51 
Gillon, M., Pont, F., Moutou, C, Bouchy, F., Courbin, F., et al. 

2006, A&A, 459, 249 
Grillmair, C. J., Burrows, A., Charbonneau, D., Armus, L., 

Stauffer, J., et al. 2008, Nature, 456, 767 
Harrington, J., Luszcz, S., Seager, S., Deming, D., & Richardson, 

L. J. 2007, Nature, 447, 691 
Horne, K. 1986, PASP, 98, 609 

Hebrard, G., & Lecavelier des Etangs, A. 2006, A&A, 445, 341 
Hubbard, W. B., Fortney, J. J., Lunine, J. I., et al. 2001, ApJ, 560, 
413 

Jolly, A., McPhate, J. B., Lecaveher des Etangs, A., et al. 1998, 

A&A, 329, 1028 
Kuchner, M. J., & Seager, S. 2005, |arXiv:astra-ph/0 504214 | 
Knutson, H. A., Charbonneau, D., Cowan, N. B., Fortney, J. J., 

Showman, A. P., et al. 2009, ApJ, 690, 822 
Knutson, H. A., Charbonneau, D., Noyes, R. W., Brown, T. M., & 

GiUiland, R. L. 2007a, ApJ, 655, 564 
Knutson, H. A., et al. 2007b, Nature, 447, 183 
Kurucz, R. L. 1979, ApJS, 40, 1 

Ku rucz R. 2006, Stellar Model an d Associated Spectra 
|http: / /kurucz. harvard.edu/grids. html) 



7 



Lecavelier des Etangs, A., Vidal-Madjar, A., Roberge, A., et al. 

2001, Nature, 412, 706 
Lecavelier des Etangs, A., Pont, F., Vidal-Madjar, A., & Sing, D. 

2008a, A&A, 481, L83 
Lecavelier des Etangs, A., Vidal-Madjar, A., Desert, J.-M., & Sing, 

D. 2008b, A&A, 485, 865 
Liang, M.-C, Parkinson, C. D., Lee, A. Y.-T., Yung, Y. L., & 

Seager, S. 2003, ApJ, 596, L247 
Liang, M.-C, Seager, S., Parkinson, C. D., Lee, A. Y.-T., & Yung, 

Y. L. 2004, ApJ, 605, L61 
Mandel, K., & Agol, E. 2002, ApJ, 580, L171 

Moutou, C, Pont, F., Bouchy, F., & Mayor, M. 2004, A&A, 424, 
LSI 

Morales-Calderon, et al. 2006, ApJ, 653, 1454 
Naylor, T. 1998, MNRAS, 296, 339 

Pont, F., Zucker, S., & Queloz, D. 2006, MNRAS, 373, 231 
Pont, F., Gilliland, R. L., Moutou, C, Cliarbonneau, D., Bouchy, 

F., et al. 2007, A&A, 476, 1347 
Pont, F., Knutson, H., Gilliland, R. L., Moutou, C, & 

Charbonneau, D. 2008, MNRAS, 385, 109 
Reach, W. T., et al. 2006, IRAC Data Handbook v3.0 
Redfield, S., Endl, M., Cochran, W. D., & Koesterke, L. 2008, ApJ, 

673, L87 

Richardson, L. J., Harrington, J., Seager, S., & Deming, D. 2006, 
ApJ, 649, 1043 



Roberge, A., Feldman, P. D., Weinberger, A. J., Deleuil, M., & 

Bouret, J.-C. 2006, Nature, 441, 724 
Seager, S.,& Sasselov, D. D. 2000, ApJ, 537, 916 
Sing, D. K., Vidal-Madjar, A., Desert, J. -M., Lecavelier des Etangs, 

A., & Ballester, G. 2008a, ApJ, 686, 658 
Sing, D. K., Vidal-Madjar, A., Lecavelier des Etangs, A., et al. 

2008b, ApJ, 686, 667 
Sharp, C. M., & Burrows, A. 2007, ApJS, 168, 140 
Showman, A. P., Fortney, J. J., Lian, Y., Marley, M. S., Freedman, 

R. S., et al. 2008, arXiv:0809.2089 
Swain, M. R., Vasisht, G., & Tinetti, G. 2008, Nature, 452, 329 
Swain, M. R., Vasisht, G., Tinetti, G., Bouwman, J., Chen, P., et 

al. 2009, ApJ, 690, L114 
Tinetti, G., Liang, M.-C, Vidal-Madjar, A., Ehrenreich, D., 

LecaveUer des Etangs, A., Yung, Y. L. 2007a, ApJ, 654, L99 
Tinetti, C, et al. 2007b, Nature, 448, 169 

Vidal-Madjar, A., Lecavelier des Etangs, A., Desert, J.-M., et al. 

2008, ApJ, 676, L57 
Vidal-Madjar, A., et al. 2004, ApJ, 604, L69 

Vidal-Madjar, A., Lecavelier des Etangs, A., Desert, J.-M., 
Ballester, G. E., Ferlet, R., Hebrard, C, & Mayor, M. 2003, 
Nature, 422, 143 

Werner, M. W. et al., 2004, ApJS, 154, 1 

Winn, J. N., et al. 2007, AJ, 133, 1828 



8 



Table 1. Limb darkening coefficients 



IRAC Channel 


Wavelength 


CI 


C2 


C3 


C4 


1 


3.6 fira 


0.6023 


-0.5110 


0.4655 


-0.1752 


2 


4.5 fim 


0.6886 


-0.8869 


0.8504 


-0.3077 


3 


5.8 fim 


0.7138 


-1.0720 


1.0515 


-0.3825 


4 


8 fim 


0.5497 


-0.8042 


0.7698 


-0.2719 



Table 2. Fitting Parameters of the Transit Curves 



Parameter 


3.6 Aim(Ch 1) 


4.5 ^jm(Ch 2) 


5.8 ^(m{Ch 3) 


8 Atm(Ch 4) 


a{Rp/FU) 


0.15449 
±0.00031 


0.15565 
±0.00033 


0.15466 
±0.00054 


0.15437 
±0.00045 


h 

a{h) 


0.6327 
±0.0099 


0.6562 
±0.0076 


0.6279 
±0.0111 


0.6651 
±0.0060 




9.136 
±0.091 


8.940 
±0.077 


9.257 
±0.107 


8.890 
±0.069 


Tq^HJD) 
ct{To){HJD) 


2454039.82881 
±0.00009 


2454428.08128 
±0.00011 


2454039.82895 
±0.00013 


2454428.08151 
±0.00009 




-0.06-0.04-0.02 0.00 0.02 0.04 -0.04-0.02 0.00 0.02 0.04 -0.04-0.02 0.00 0.02 0.04-0.04-0.02 0.00 0.02 0.04 
Phase 



Fig. 1. — Transit light curves, fits and residuals for each channel. Data are not rebinned. Typical Icr error bars on individual measurements 
are represented. The raw weighted light curves (top panel) have to be corrected for large fluctuations correlated to the 'pixel phase' and 
to the ramp baseline. The corrected light curve is plotted in the middel panel without the rejected points (rejection at more than 3cr ). 
Overplotted is the fit with limb darkening taken into account. The bottom plot show the residuals from the best fits. 
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Fig. 2. — Channel 4. Top panel: Rp/Ri, ratios extracted when using linear, polynomial of a 3rd degree, and logarithmic baselines as 
function of the number of first photometric points removed from which the transit light curve is fitted. In the abscisse, point removed 
correspond to the full transit light curve (TLC). The left over TLC when padded at a phase of —0.03 correspond to the removal of ~ 430 
measured points. The diamond point on the right hand side with its error bar, indicates the mean value with its error bar obtained by 
bootstrap when using a logarithmic baseline (Sect. I3.5|l . Bottom panel: Root-mean-square (RMS) resulting from the fit using the three 
different baselines describe above as function of first exposures removed. The fit is always better when u sing the logarithmic function. This 
is particularly true at the beginning of the time series, where the ramp effect is the strongest (See Sect. |3!T1 . 
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Fig. 3. — Root-mean-squarc (RMS) of binned residuals versus bin size for ehannel 1 for two different pixel-phase effect decorrelation. 
The solid line is proportional to N^^^^. Top panel: RMS resulting from the fit with f parameter used to fit the corrective coefficient in 
the X and Y direction axis. Bottom panel: RMS resulting from the a quadratic fit with 5 parameters. Negligible systematics appear in the 
residuals when using 5 parameters to characterize the pixel-phase effect. 
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Fig. 4. — Measured values of Rp/Ri, ratios in the four Spitzer /IKAC channels for HD 189733b. Results in channel 1 at 3.6 /jm and in 
channel 3 at 5.8 ^tm are from the same data sets but from different analysis (Ehrenreich et al. 2007, Beaulieu et al. 2008, and this work). 
Our analysis presented in this paper is more robust against systematics than the previous ones. New observations in channel 2 at 4.5 fira 
and in channel 4 at 8.0 /xm are also presented here. Results from our present analysis are plotted with red squares. Previous estimates 
from the same dataset at 3.6 ^m and at 5.8 /^m are plotted with green dots (Ehrenreich et al. 2007) and with blue triangles (Beaulieu et 
al. 2008). An independent measurement from Knutson et al. (2007a) at 8 ^tm in Channel 4 is plotted with a purple star. The Rp/Rt ratio 
obtained at 4.5 fim is larger (4(t) than expected with either water absorption only (light blue line) or with Rayleigh scattering by small 
particles (grey dotted line; Pont et al. 2008; Lecavelier des Etangs et al. 2008a). Therefore, other species must be present in the atmosphere 
of the planet absorbing at this wavelength. This supplementary absorption could be due to the presence of CO molecules in the planetary 
atmosphere. 



